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We extract all the invariants (i.e. all the functions which do not depend on the choice of phase- 
space coordinates) of the dynamics of two point- masses, at the third post-Newtonian (3PN) ap- 
proximation of general relativity. We start by showing how a contact transformation can be used 
to reduce the 3PN higher-order Hamiltonian derived by Jaranowski and Schafer Ei to an ordinary 
Hamiltonian. The dynamical invariants for general orbits (considered in the center-of-mass frame) 
are then extracted by computing the radial action variable <j> p r dr as a function of energy and angular 
momentum. The important case of circular orbits is given special consideration. We discuss in detail 
the plausible ranges of values of the two quantities astatic, ^kinetic which parametrize the existence 
of ambiguities in the regularization of some of the divergent integrals making up the Hamiltonian. 
The physical applications of the invariant functions derived here (e.g. to the determination of the 
location of the last stable circular orbit) are left to subsequent work. 



I. INTRODUCTION 



Binary sytems made of compact objects (neutron stars or black holes) are the most promising candidate sources 
for ground based interferometric gravitational- wave detectors such as LIGO and VIRGO. In the case of (stellar) black 
holes the gravitational waveform will enter the detector bandwidth only at the last stage of the inspiral motion, 
just before the inspiral turns into a plunge. The detection of these gravitational signals will be possible only by 
correlating the detector output with sufficiently accurate copies of the expected signals ("matched filtering"). It is 
therefore very important to have the best possible analytical control of the dynamics of general relativistic two-body 
systems. The equations of motion of binary systems have been derived some years ago at the 5/2 post-Newtonian 
(2.5PN) approximation^ Q. Recently, it has been possible to derive the third post-Newtonian (3PN) dynamics of 
point-mass systems [||, though with some remaining ambiguity due to the need to regularize the divergent integrals 
caused by the use of Dirac-delta-function sources. The knowledge of the conservative part of the dynamics of binary 
systems must then be completed by a correspondingly accurate knowledge of the gravitational-wave luminosity (used, 
heuristically, to derive an accurate estimate of the radiation damping effects which drive the inspiral motion of binary 
systems). The gravitational-wave luminosity is currently known to the fractional 2.5PN accuracy ||. By combining 
the 2PN-level conservative dynamics |l| with the 2.5PN gravitational- wave luminosity || one has recently constructed 
some (improved) filters (P-approximants) for application to gravitational- wave data analysis problems. 

The main purpose of this work is to extract all the invariants, i.e., the functions which do not depend on the choice 
of coordinates (in space or phase-space), of the 3PN dynamics derived in Ref. 0. This task is important for three 
reasons: (i) some of the invariant functions are directly useful for deriving the 3PN "phasing formula" of inspiralling 
binaries, i.e., for constructing 3PN-accurate gravitational- wave filters; (ii) we shall use, in a companion paper [0, 
some of the invariants derived below to determine the location of the Last (circular) Stable Orbit which marks the 
transition between the inspiral and the plunge, and (iii) from a practical point of view, the dynamical invariants will 
be quite useful for comparing the 3PN ADM Hamiltonian dynamics of [pj with forthcoming derivations of the 3PN 
equations of motion in harmonic coordinate systems B . 



1 We recall that the "nPN approximation" means the obtention of the terms of order (v/c) 2n ~ (Gm/c 2 r) n in the equations 
of motion. 
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II. REDUCTION OF THE 3PN HIGHER-ORDER HAMILTONIAN 



It was shown some years ago that, in most coordinate systems, the conservative part of the PN-expanded 
equations of motion of two body systems, say x a = A. a (x b , x ), where a, b — 1,2 and where A. — Ao + c~ 2 A2 + 
c _4 A4 + c~ 6 Aq + . . . , do not follow from any ordinary Lagrangian L(x a , x a ). For instance, in harmonic coordinates, 
and at the 2PN level, one needs to consider an acceleration-dependent Lagrangian L2p™ omc (x a , x a , x a ) [pj. However, it 
was shown in Ref . , and, more generally, in Ref . JtJ , that any higher-order PN-expanded Lagrangian L (x , x a , x a , . . . ) 
(where higher derivatives enter only perturbatively) can be reduced to an ordinary Lagrangian L'(x' a , ~k! a ) by a suitable 
(higher-order) contact transformation x' a (t) = x a (t) — e a (x ,Xf,, . . .). At the 2PN level the class of coordinates where 
the dynamics admit an ordinary Lagrangian is rather restricted |?]] , but it includes in particular the ADM coordinates 
P,^0|. At the 3PN level, Jaranowski and Schafer, who worked within the ADM canonical formalism, have found 
that the (conservative) ADM dynamics could not be derived from an ordinary Hamiltonian H(x a ,p a ) (equivalent 
to an ordinary Lagrangian L(x a ,x a )), but that instead it could be derived from a certain higher-order Hamiltonian 
H(x a , p a , x a , p a ). [This higher-order matter Hamiltonian is defined by eliminating the field variables hj^ , hj? in a 
certain Routh (i.e., mixed) functional i?(x a , p a , hj?, hjj 1 )', see Eq. (33) of Ref. Q.] The meaning of this higher-order 
matter Hamiltonian H(x a , p Q , x a , p a ) is that the equations of motion of the matter can be written (after elimination 
of the TT variables) as 

. _ SH(x b ,p b ,x b ,p b ) 5H(x b ,p b ,x b ,p b ) 

x a — r i Pa — r i V^-JJ 

OPa 0X Q 

where 6/5p a and 6/6x a denote functional derivatives: 

(2.2) 





It is easily seen that the Hamilton-like equations of motion (2.1) are equivalent to the Euler-Lagrange equations 
derived by extremizing the action functional 

S[x a (i),p a (t)] = J diZ(x b ,pb,x b ,pb), (2.3) 

where 

L{x b , p b , x b , p b ) = p a x a - H(x b , p b , x 6 , p b ). (2.4) 

Indeed, we have 

SL SH 5L . SH rn . 

— = x a - — , — = -p a - — . (2.5) 

OPa OPa 0X a dX a 

To derive the dynamical invariants of the 3PN equations of motion it is convenient to introduce new coordinates 
(in phase-space (x a ,p a )), 

x' a =x a + 5x a , p' a = Pa + 5p a , (2.6) 

such that the equations of motion for (x^pJJ become some ordinary Hamilton equations 

_ dH'{x' bl p' b ) dH'(x' b ,p' b ) 
a ~ dp' a ' dx' a ■ {2 - 7) 

This is equivalent to requiring that the (x„, p^) equations of motion derive from an action functional of the form 

S'[x' a (t),p' a (t)] = J dtL'{xlp' b ,xlp' b ), (2.8) 

where^ 



2 As defined L' does not depend on p', but, by symmetry, it is convenient to allow for such a dependence (which can be easily 
introduced by transforming p'x' by parts). 
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L'(x;, P ;,x;,pa = P x-ff'(4,pl). 



(2.9) 



When so expressed at the level of action functionals, the problem of reducing the 'higher-order' action S[x a (t), p a (t)] 
to an ordinary action S"[x' a (t), p a (£)] is quite similar to the problem of the order-reduction of (perturbative) higher- 
order actions S^x^i)] = J dt L(x a ,x a ,x a , . . .) which was solved in full generality in Ref. [Q. It is then an easy task 
to adapt the techniques used in Q| to solve our present problem. 

First, we note that, taking into account the fact that the Hamiltonian has the perturbative structure, 



H(x a ,p a ,± a ,p a ) = H (x a ,pa) + c 2 iJ 2 (x a ,p a ) + c 4 i/ 4 (x Q ,p a ) +c 6 H 6 (x Q , p a , x a , p Q ) , the identities (gj) for 
the variational derivatives of L(x a , p a , x a , p a ) imply 

SL SL 
x a = v 2PNa (x, p) + — + 0{c- 6 ) = v Na (x, p) + — + 0(c- 2 ), (2.10a) 

OPa 0p a 

SL SL 
Pa = q 2 PNa(x,p) - — +0{c- 6 ) = q N a(x,p) - — + 0{cr 2 ). (2.10b) 
ox Q 0x a 

Here V; a and q^ a (with i — 2PN or N) denote some explicit functions of x and p which are the right-hand-sides of the 
usual Hamilton equations at the 2PN (or Newtonian) level. For i nstan ce, at the Newtonian level vn q (x, p) = p /m , 
qNa(x, p) = — Gm a 77ib(x a — x )/|x a — x | 3 . Note that in Eqs. ( 2. 10|) , and in the reasoning below, we are working 



'off shell', i.e., we consider virtual motions which do not necessarily satisfy the equations of motion. [As we are 
using identities for action functionals, it is essential to work off shell.] Inserting the perturbative identities ( [2.1C ) in 
H (x, p, x, p) and Taylor expanding yields the identity^ 

QJJ fij^ QL[ SL 

#(x,p,x,p) = i?(x,p) + — -7— — - — (- double-zeros +C(c" 8 ), (2.11) 

ax op op ox 

where H(x, p) denotes the naive 'order-reduced' Hamiltonian obtained by using the (lower-order) equations of motion 
to eliminate the higher-order derivative terms (a wrong procedure in general): 

#(x, p) = #(x, p, v N (x, p), q N (x, p)) + 0(c- 8 ). (2.12) 

[As indicated in Eq. ( [2.12] ), it is sufficient to u se the Newtonian equations of motion because x and p enter only at 
the c -6 level.] The 'double-zeros' in Eq. ( [2.11 ) denote all the terms generated by the Taylor expansion which would 



be at least quadratic in SL/Sx and SL/Sp. As is well known [|ll],[7| such terms can be perturbatively neglected (even 
off shell) because they do not contribute to the equations of motion at the 3PN level. 



If one inserts the identity (2.11) in the original action functional (2.3), one sees that 5[x, p] has the desired ordinary 
form ( |2.8| ) modulo some extra terms which are linear in 6L/5x and SL/Sp. As in Ref. such terms can be eliminated 
by noticing that (to first order) the effect of the shift of dynamical variables ( |2.6| ) on the Lagrangian L(x, p,x, p) 
reads (by virtue of the definition of functional derivatives) 

SL SL d 

Z(x, p, x, p) = Z(x', p', x', p') - — Sx - —Sp + -Q(Sx, Sp), (2.13) 

Ox op dt 

where Q(o~x, Sp) is some linear form in <5x and 0~p. [It is sufficient to work to linear order because <5x and Sp are 
0(c" 6 ).] 



By combining the two identities (2.11) and (2.13) we find that, if we define the ordinary Lagrangian (considered in 



phase-space) L, associated to the ordinary (naively reduced) Hamiltonian H(x, p), Eq. (2.12), 

L(x,p,x,p) = p Q x a - i?(x,p), (2.14) 

we have the identity 



In Eq. (2.11) and elsewhere, we abbreviate the notation by suppressing the (summed over) indices in x = x l a , p = p a i, etc. 
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7/ • -n t, / / •/ • /\ SL (. 8h\ 6l{. 8h\ 
L(x,p,x,p)=L(x,p,x,p)--l*x- — l-_fip+_ I 



— Q(5x, 5p) + double-zeros +C(c 8 ). 



Therefore, if we shift the phase-space coordinates by defining 

£>c„ = — . p a - p a 



dH 

dpa'' 



dH 
3x a ' 



(2.15) 



(2.16) 



we find (noticing tha t bo th total differentials, and double-zeros, are negligible in action functionals) that the original 
equations of motion (2.1) are transformed, when rewritten in (x',p') coordinates, into the Euler-Lagrange equations 
of the 'ordinary' phase-space Lagrangian L(x', p', x', p'), i.e. into ordinary Hamilton equations 



Off(x',p') 
dp 1 



+ 0(c~ 8 ), 



aff(x',p') 



+ 0(c~ s ). 



(2.17) 



Summarizing: it is licit to naively reduce the higher-order original Hamiltonian ff (x, p, x, p) by replacing the lower- 
order equations of motion in the offending derivatives x, p (thereby defining the reduced Hamiltonian ff(x, p), Eq. 
( [2.12 )), if one adds the correcting information that the ordinary Ham ilton equations defined by the reduced if (x, p) 
hold in the new phase-space coordinates (x',p') defined by Eq. ( 2.16| ). 

Let us note for completness that the conserved energy £' = if(x', p') defined by the autonomous Hamiltonian H 
becomes, when expressed in the original variables 



where 



& = if(x', p') = if(x, p) + - = £ + O(o-*) 

ox op op ox 



~, ... .dH .dH 
if (x, p, x, p - x— - p— 
ax op 



(2.18) 



(2.19) 



is, indeed, easily checked to be a conserved quantity associated with the higher-order dynamics (2.1). 

Let us now consider the explicit application of the previous results to the case at hand. Following [|2 it is sufficient to 
consider the dynamics of the relative motion of a two-body system, considered in the center-of-mass frame (pi + P2 = 
0). It is convenient to work with the following reduced variables^ 



xi - x 2 
GM : 



= El El 



GM' 



if 



XR 



ff NR 



where 



M — mi + 7712, (i 



m\rri2 
M 



M (mi + ?7i2) 2 



(2.20) 



(2.21) 



In Eq. (2.20) the superscript NR denotes a 'non-relativistic' (higher-order) Hamiltonian, i.e. the Hamiltonian without 

R is, to sta 

dr d(xi — x 2 ) 



the rest- mass contribution Mc 2 . ff NR is, to start with, a function of r, p, r, and p. Here the dot denotes the 
reduced-time derivative: 



dt 



dt 



dp G dpi G dp2 
^ dt v dt v dt 



(2.22) 



From Eq. ( fy.lt ) above and Eq. (68) of J || one finds that the coordinate shift needed to transform if NR (r, p, r, p) into 
an ordinary Hamiltonian reads (on shelly, i.e. by using the (Newtonian) equations of motion after the differentiations 
exhibited in Eq. fl2.16|)) 



4 Note that [g] use units where 16nG = 1. 

5 Working on shell here is equivalent to neglecting some double- zero terms, see, e.g. Ref. 
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6r = r' — r = 



g H NR 



dp 



on shell 



^ 2 ((n • p) 4 - (p 2 ) 2 ) + i((19 - 26u)u (n • p) 2 - 5(13 + 4z>p 2 )l 



n 



1 

+ 12 
5p = p' p 



^ 2 (3p 2 + (n • p) 2 )(n • p) - (14 - v)v (n • p)- 



P , 



(2.23a) 



dr 



on shell 



^ 3 ((p 2 ) 2 + 2(n • p) 2 p 2 + 5(n • p) 4 )(n • p)l 



-(3(1 - 9is)vp 2 + 5(17 - 2<dv)v (n • p) 2 )(n • p)4 + 7T^( 11408 - 9157r2 + 36 ( 48 - ^)v)v (n • p)^ 



48 



n 



^ 3 ((p 2 ) 2 + 2(n • p) 2 p 2 + 5(n • p) 4 )^ + ^((4 - 15^p 2 + (6 + lv)v (n • p) 2 )-^ 



1536 



(15(61tt 2 



4(45tt 2 - 592)i^- 



(2.23b) 



Note again that 8r and 5p are of 3PN order. 



III. ORDER-REDUCED 3PN HAMILTONIAN 



As explained above the phase-space coordinate transformation (2.22) maps the original higher-order 3PN ADM 



dynamics onto ordinary Hamilton equations, with an Hamiltonian H (r', p') defined by the 'naive' procedure (2.12), i.e. 

by using in i? NR (r, p, r, p) the following replacement rules (in our reduced units where -ffNewtonian( r ' P) = P 2 /2 — 1/t, 
with r = Irl) 



r^r, p^p, r^p, p 



r 

1/3 ' 



(3.1) 



To simplify the writing we shall henceforth drop all the primes (but the reader should remember that we henceforth 
work in the shifted coordinate system (r', p')). We also drop the overbar, on the Hamiltonian, indicating that we have 
used the replacement rules (3J). Finally, using Eqs. (68) and (71) of we obtain the following ordinary Hamiltonian 
(with n = r/r) 



H NR (r, p) = H N (r, p) + ^H 1PN (r, p) + ^H 2 pn (r, p) + ^H 3 pn (r, p) 



(3.2) 



where 



Hn (r,p) = , 

2 r 

iT 1PN (r, p) = i(3«/ - l)(p 2 ) 2 - \ [(3 + „)p 2 + u(n • p) 2 ] \ + ± 



(3.3) 
(3.4) 



i?2PN (r, p) = 1 (1 _ 51/ + 5^ 2 ) (p 2 ) 3 + i [(5 - 20u - 3^ 2 ) (p 2 ) 2 - 2^ 2 (n ■ p) 2 p 2 - 3*/ 2 (n ■ p) 4 ] 1 



i[(5 + 8,)p 2 + 3Kn- P ) 2 ]l-i(l + 3,)l 



(3.5) 



i? 3 PN (r, p) = 4s (-5 + 35*/ - 70j/ 2 + 35t/ 3 ) (p 2 ) 4 



+ 77 [(- 7 + ^ - 53 ^ 2 - 5 ^ 3 ) (P 2 ) 3 + ( 2 - 3*V(n • p) 2 (p 2 ) 2 + 3(1 - jV(n • p) 4 p 2 - 5z/ 3 (n • p) 6 ] - 
16 r 
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1 

16 



-27+ 136i/ + 109^ 



(p 2 ) 2 + ^(17 + 30^(n.p) 2 p 2 



1 

12 



(5 + 43i/)j/(n • p) 4 



25 . 1 2 335, 55 , 

h ( 7T )V V' 

8 ^64 48 ' 12 



85 
16 



64 



27 



v(n-p) : 



1 .109 21 2 , 

- + ( ir 2 )v 

8 v 12 32 ' 



- ^kinetic (p 2 - 3(n ■ p) 5 



' ^static" 



(3.6) 



The parameters Wkinetic and w s tatic appearing in the 3PN Hamiltonian i?3PN parametrize the existence of ambiguities 
in the regularization (as it is presently performed) of some of the divergent integrals making up the Hamiltonian. 

It was shown by Damour |12| that the internal structure of the compact bodies making up a binary system (e.g. 
a neutron star versus a black hole) start influencing the dynamics only at the 5PN level (G 6 /c 10 ; see Eq. (19) in 
Sec. 5 of Ref. @). Therefore, one expects that the present 3PN ambiguities are of technical nature, and not linked 
to real physical ambiguities. The best way to resolve these ambiguities would be to work out in detail, at the 3PN 
level, the general matching method outlined in Ref. jl2), and applied there at the 2PN level. However, the technical 
difficultities of working at the 3PN level are such that this route has not yet been attempted. 

The 'kinetic ambiguity' widnctic^ 2 (p 2 — 3(n • p) 2 ) /r 3 (proportional to the kinetic energy p 2 ) was explicitly intro- 
duced in Ref. g] (see Eq. (69) there; with cj thcrc corresponding to uj^^ iir modulo the shift of normalization discussed 
below). The 'static ambiguity' w s tatic*V^ 4 was discussed in Ref. [ fl3| . In the present paper, we normalize the static 
ambiguity such that the value w s t a tic = corresponds to the static interaction^ between two 'Brill-Lindquist black 
holes' (considered instantaneously at rest; see |l3|]). The normalization of the kinetic ambiguity used in Eq. (3.6) 
of this paper also differs from that used in Eq. (71) of The motivation for this new normalization is given in 
Appendix A of the present paper, where we calculate, following the route of Ref. a new 'reference' Hamiltonian 
which corresponds to this normalization. In terms of these two normalizations, the Hamiltonian explicitly written 
in |J2) corresponds to the values astatic — 1/8 and Wkinctic = 7/16 + u>] this means that (after having applied the 
replacement rule (3.1) and dropped the primes) 



H. 



this paper 
3PN 



(r,p) 



#3PN ( r > P) 



1 



^static 
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^kinetic — yg — ^ 



(p 2 -3(n.p) 5 



(3.7) 



Though the two parameters w s t a tic and osmetic appear separately in Eq. (B.6), all the dynamical invariants of the 
3PN Hamiltonian involve only the combination 



a(v) = ^W s tatic + ^ ^kinetic- 



(3.8) 



Our results below will verify explicitly that only the combination a appears, but it is, in fact, easy to understand why. 
Indeed, as we shall calculate the dynamical invariants of H(r,p), our results would be the same if we transformed 
H by an arbitrary canonical transformation. If we consider 'a small' canonical transformation of order 1/c 6 , it is 
enough to work at linear order. At linear order a canonical transformation is parametrized by an arbitrary generating 
function g(r, p), and its effect on the Hamiltonian is to change it by the Poisson bracket {g, H}: 



SH(r,p) = {g,H} 
If we choose g(r, p) = Ar • p/r 3 we find 



dg dH dg dH dg x l dg 



dx l dpi dpi dx l 



' dx l 



dpi 



0(g/c 2 ). 



SH 



(p 2 -3(n.p) 2 



^4' 



(3.9) 



(3.10) 



Therefore by choosing (for instance) A = — ^ 2 Wkinotic, we see that we can eliminate osmetic at the cost of replacing 

^static by ^Wstatic ~ A = I^ s tatic + V Wkinetic- 



An important issue, when discussing the consequences of the 3PN Hamiltonian (3.2), is the influence of the lack 
of precise knowledge of the combination cr, Eq. (3.8), on physical observables. First, let us emphasize that, among 



By 'static interaction' we mean the part of the Hamiltonian which remains when setting to zero both the momenta pi, P2 
and the independent gravitational- wave degrees of freedom hj^ , 7r l:,TT (time-symmetric conformally flat data) . 
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the hundreds of contributions to the 3PN Hamiltonian which have been computed in Ref. Eq|, most of them, though 
they are in general given by a divergent integral, seem to be regularizable in an unambigous manner. Indeed, if one 
considers separately the divergences near one particle (say in a volume V\ near particle 1), and regularizes them 'a la 
Riesz', i.e. by introducing a regularizing factor (|x — Xi|/Zi) ei = (ri/Zi) Cl in the integrand, the analytic continuation 
in £1 of the regularized integral ii(ei) = J v d 3 x (ri/Zi) ei F(x, xi, X2) does not contain a pole oc 1/ei for most terms. 
This suggests that the wcll-dcfincd analytic continuation of Ji(ei) down to t\ = defines, without apparent ambiguity, 
the regularized value of J v d 3 xF(x, xi,x 2 ). [In the notation of Appendix B of § we are considering that v JS = 

(i.e. £2 = ^ JS e = 0), in a Vi-restricted integral which behaves as Ai + Ci^ JS //i JS = A\ + C^/ei, before taking the 
limit ei = ^ Js e — * 0.] 

However, a limited class of 'dangereous' integrals involve a simple pole as e% — > 0: ii(ei) = Zi^e^ 1 + ln(i?i/Zi)) + A\ 
(where the 'infrared' lengthscale R\ is linked to the size of the volume Vi). A remarkable fact (emphasized in Section 
IV of ||) is, however, that the complete combination of dangerous integrals appearing in i?3PN is such that all the 
pole terms cancel: = {). This also implies that the logarithms depending on the 'ultraviolet' regularizing length 

l\ cancel in the Hamiltonian In other words, the combination of dangerous integrals appearing in i?3PN is either 
'finite' (convergent) or globally of a 'non dangerous' type. This is good news, but this still leaves an ambiguity in the 
finite value of ^ because it was noticed in [0 that the final regularized result depends on the way one writes 

the 'bare' integrand -F(x, xi, X2) (when transforming it using standard rules: integration by parts, and Leibniz's rule). 
In addition to this ambiguity linked to a subset of divergent integrals, there are also ambiguities in the 'contact' 
terms in the Hamiltonian (of the type J d 3 x8(x — XjJS^x, xi, X2)). One needs to regularize the (in general singular) 
limit lim x ^ Xl S*(x, Xi,X2). In [^) one used an Hadamard-like 'partie-finie' prescription (Pfi) for doing that. This 
prescription is, however, ambiguous at 3PN, notably because Pfi(4>f 2 )) [Pfi(0(2))] 4 where 0( 2 ) denotes the 
Newtonian potential. As said above, the final amb igui ties in the Hamiltonian concern only two quite specific types of 



terms, parametrized by w sta tic and akinetic in Eq. (3.6) 



In Appendix A of this paper, we discuss in more detail these ambiguities, and, we try to estimate what are the 
corresponding plausible ranges of values of the two ambiguity parameters ^static and ^kinetic- For instance, when 
transforming (e.g. by integration by parts) the divergent integrands we generate some (rational) values for w s tatic and 
^kinetic- By doing this in various ways, we get an idea of the range of values that these parameters are likely to be in. 
The final result of Appendix A below is 

- 10 < Static < 10, -10 < W kino tic < 10. (3.11) 



As the symmetric mass ratio v = 17111712/ (mi + m2) 2 ranges between and 1/4, the ranges (3.11) imply that —3.125 < 
o(y) < 3.125. In the following, we shall therefore consider as fiducial range of variation for the (not yet known) 
combination a the simple range 

- 3 < -12^ < o(v) < Ylv < 3. (3.12) 



IV. DYNAMICAL INVARIANTS OF THE 3PN HAMILTONIAN 



In this section we derive a complete (and, in fact, overcomplete) set of invariant functions of the 3PN dynamics. 
Only such functions enter the principal observables that can be measured from infinity (say, by gravitational wave 
observations). However, not all the invariants we derive here can be considered as being directly observable. In a 
subsequent paper, we shall discuss in detail one of the most important observable: the Last Stable (circular) Orbit. 

We follow closely a work of Damour and Schafer jli} in which they derived the invariants of the 2PN dynamics. One 
first uses the invariance of the Hamiltonian under time translations and spatial rotations. This yields the conserved 
quantities 

cNR t 

E^ — ^# NR (r,p), jE — =rxp. (4.1) 

Here £ NR denotes the total 'non relativistic' energy (without rest-mass contribution), and J the total angular mo- 
mentum of the binary system in the center of mass frame. Using the Hamilton- Jacobi approach, the motion in the 
plane of the relative trajectory is obtained, in polar coordinates (r, <j)), by separating the variables t = tj (GM) and </> 
in the reduced planar action 

S 



S = = -Et + j<j> + / dry/R(r,E,j). (4.2) 
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Here j = Jjl and R(r, E, j) is an 'effective potential' for the radial motion which is obtained by solving the first 
equation (hi) for p 2 = (n • p) 2 after having replaced p 2 by 



p 2 = (n.p) 2 + (nxp) 2 =p 2 + ^. (4.3) 

Working iteratively in the small parameter c~ 2 , and consistently neglecting all terms 0(c~ 8 ), one finds that R(r, E, j) 
is given, at 3PN order, by the following seventh-degree polynomial in 1/r: 

2B C nlPN D 2PN D 2PN D 3PN D 3PN 

i?(r , jB)i)=A+ _ + _ + ^_ + ^_ + ^_ + ^_ + ^_. (4 .4) 

The coefficients A, B, C start at Newtonian order, while the extra terms Di/r l+2 start at the PN order indicated 
as superscript. All the coefficients A, B, C, Di are polynomials in E and j 2 . Their explicit expressions are given in 
Appendix B. 

The Hamilton- Jacobi theory shows that the observables of the 3PN motion are deducible from the (reduced) 'radial 
action integral' 



i r (E,j) = — dry/R{r,E,j). (4.5) 

For instance, the periastron-to-periastron period P and the periastron advance A<i> = $ — 2tt per orbit are obtained 
by differentiating the radial action integral: 

P d 

ME, 3), (4.6) 



2itGM dE 
$ A$ d 

- = 11- = --^,! (4.7) 
2ir 2ir oj 

To evaluate the invariant function i r (E,j) we follow Ref. which was itself following a method devised by 

Sommerfe ld jllq l within the context of the 'old theory of quanta'. This method is explained in detail in the Appendix 
B of Ref. |l4|]~The basic idea is the following: define Q(r) = A + 2B/r + C/r 2 and consider the integral 



Here e is a formal expansion parameter (actually in the final calculations, one takes into account the fact that the 
higher Z?„'s are multiplied by e 2 or e 3 , with e ~ 1/c 2 being in factor of Z?} PN ) and one wishes to compute the expansion 
of 7(e) as e — > 0: 1(e) = Iq + e/i + e 2 h + +0(e 4 ). To start with, the limits of integration, n (e) and r 2 (e), are the 
two exact (e-dependent) real roots of R(r) — Q(r) + e^2 D n /r n+2 . [We work in the case where A < 0, B > 0, and 
C < 0.] The idea is to consider the function \J R(r) as a complex analytic function defined in a suitably cut complex 
r-plane: with, say, a cut along the real segment lin king r\(e) to ^(e), and additional cuts from, say, — oo to the other 



roots. This allows one to rewrite the real integral (4.S) as a contour integral 



1 



1(e) = —& dry/Rj^e), (4.9) 
^ Jc 

where one should note the different numerical prefactor, and where C denotes any (counterclockwise, simply winding) 
contour enclosing the two real roots r±(e), T 2 (e), and avoiding all the cuts. [The phases of r — r% and r — r 2 are defined 
to be zero when r lies on the real axis on the right of r 2 , while the phase of \f~A is defined as +tt/2.] If we keep the 
contour C a way from the segment [ri(e), r 2 (e)], we can now directly expand the integrand yj R(r, e) in the contour 



integral (4J) in powers of e. This generates the successive terms of the expansion 1(e) = Iq + el\ + e I 2 + ■ ■ ., with 



each term /„ being given by a contour integral made of a sum of contributions of the form 

If,/. 2B C\^ 9 I 



J^ = ^f c *r{A + - + ^) -. (4.10) 

Each basic integral J q , m appears in the expansion of / multiplied by a coefficient which is a polynomial in the D n 's. 
[At any given PN order, there are only a finite number of integrals J q ^ n to compute; see below.] 
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There are then two ways to compute the J^m's. The simplest (the one advocated by Sommerfeld) is to expand out 
the contour C (away from the natural cut [n(0), 7*2(0)] associated to the 'unperturbed' quadratic form Ar 2 + 2Br + C) 
until it is deformed into: (i) a clockwise contour Co around the origin r = 0, and (ii) anticlockwise contour Coo at 
infinity. [Here, one considers that the J-integrals are defined in a newly cut plane, with a cut only along the segment 
[t*i(0), 7*2(0)].] The value of J 9jm is then simply given by applying (twice) Cauchy's resi due t heorem: it is enough to 



read off the coefficients of \/r in the two Laurent expansions of the integrand of Eq. ( 4.10 ) for r — * and r — ► 00 
(keeping track of phases and signs). The other way to compute the J's consists in shrinking down the contour C onto 
the real axis, so as to get (twice) a real integral from 7*1 (0) + 77 to 7*2(0) — 77 (with 77 > 0), plus two circular integrals 
around ?*i(0) and 7*2(0). As shown in Appendix B of |14| one can then prove that the e-expansion (Exp e ) of /(e) is 
simply given by 



Exp £ (/(e)) = Pf 



2 

2~/7 



r 2 (0) 



i-i(O) 



drExp e (yR(r, e)) 



(4.11) 



where Pf denotes Hadamard's partie finie of the real r-integral. Each integral appearing in the formal e-expansion 
of \/R{r, e) on the right-hand side of (4.11) is again the same combination of J-integrals as above, but now the J- 
integrals are real integrals along the segment [?*i(0), 7*2(0)], of which one must extract their 'partie finic'. The parties 
furies of all those real integrals are easily computed, e.g., by using any simple trigonometric paramctrization of the 
radial variable (r = <z(l — e 2 )/(l + e cost;) orr = a(l-e cosu)) to compute the indefinite integrals J dr[Q(r)] 2 -i r ~ m _ 

The only integrals that one needs to compute at 3PN order are J q ^ m with < q < 3: when q = 0, 771 = 0; when 
q = 1, 3 < m < 7; when q = 2, 6 < m < 8; and when q = 3, m = 9. We have computed these integrals by the two 
methods indicated above, and checked that the results agree. The results are given in Appendix B below. 

Finally, one must reexpand all the coefficients A, B, C, D n (which contain various orders in 1/c 2 ) to get the 
complete 3PN expansion of the radial action i r (E,j). The final result reads 



13 1 

*r( E >j) = -j{ 1 -#p-j 



1(35- 



. 1 1, ,E' 
10i/)^j + 5(15-61/)^ 
3 2 f 



3 1 E 

-{i x {u) - <t(is))- + (i 2 (u) - a(^)- + M 3 (u) — 

A J .1 J 



1 



1 1 



1 1 



1 



1 + ^(15 - „) E + -,-(35 + 30, + 3, 2 ) - & 



(4.12) 



where 



ii(u) 



77 



41 
— 7 

64 



125 



83 



, v H v 

3 24 



(4.13a) 



ii(y) 



105 / 41 9 218 

7T I V 

64 I 



h(v) - i (5-5i/ + 4i, 2 ), 



221 
; 

24 



128 



(21-105i/+15, 2 + 5, 3 ). 



(4.13b) 

(4.13c) 
(4.13d) 



The differentiation (Elj) then leads to the following expression for the periastron-to-periastron period P 



1 



2ttGM (-2£) 3 / 2 [ c 2 4 



h 5 _ 2v ) { 2E)V2 - A(35 + 30, + 3, 2 ) E 2 
2 j 32 



(i 2 (i/) - o-(v)) { 2 f 3/2 - 3z 3 (,) ( 2 ^ )5/2 + 5*4(1/) E 3 



(4.14) 
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wher e i^iy), i%{v), and ia{v) are given in Eqs. (4.13b), (4.13c), and (4.13d), respectively. Similarly, the differentiation 
(4.7) yields, for the dimensionless parameter 



"2T 



$ - 2tt 
2vr 



measuring the fractional periastron advance per orbit, the result 



k = — — <1 + — 

c 2 j 2 \ c 2 



~(7-2*/)i + i(5-2*/)£ 



1 



§(i x (i/) - <t{v))± + (i 2 {v) a(u))^ + 3i 3 (u) E 2 



(4.15) 



(4.16) 



where ii(v), i2{v), and £3(1/) are given in Eqs. ( 4.13a ), ( 4.13b| ), and ( 4.13c ), respectively. 

We shall also follow Ref. Wm in giving the Hamiltonian as a function of the Delaunay (reduced) action variables 



■ 3 = 



fiGM' 



HGM ' 



m = j z 



fiGM 



(4.17) 



In the quantum language, JV/h is the principal quantum number, J/h the total angular-momentum quantum number, 
and J z /Ti the magnetic quantum number. [They are adiabatic invariants of the dynamics and they are (approximately) 
quantized in integers.] By rot ation al symmetry, the (reduced) magnetic quantum number m does not enter the 
Hamiltonian. By inverting Eq. ( 4. IS ) one finds 



H{n,j,m) = -r^< 1 + 4 



jn 4 n z 



|( 7 - 2 ^)i 
2 j°n 



■^L _ | (35 - 4«/)^ + 1(145 - 151/ + v 2 )-. 



3( <lM - .(,))^ + §(7-2,)^ 



(^M-^))^3-f (20-Si/)^ 



1 



1 



-(275 - 50j/ + 4^)— - —(6363 - 805i/ + 90i/ 2 - — 



64 



(4.18) 



where 



303 



41 
64' 



1427 
12 



281 
/ 

24 



(4.19) 



The results (|4.14|) , ( ^jj ), and (J4JJ) confirm the 2PN results of |14j and extend them to the 3PN level. 

The angular frequencies of the (generic) rosette motion of the binary system are then given by differentiating H 
with respect to the action variables. Namely, 



2n 



Radial - p - QM 



1 dH(n,j,m) 



dn 



(4.20) 



periastron 



2nk 



1 dH(n,j,m) 



GM 



dj 



(4.21) 



Here, w ra diai is the angular frequency of the radial motion, i.e., the angular frequency of the return to the periastron, 
while Wperiastron is the average angular frequency with which the major axis advances in space. 

It is interesting to note that the ambiguity par amet er a{y) enters only in two of the six independent com bina tions 
of n and j which enter at the 3PN level, cf. Eq. ( 4.1 8| ) . We note also that if we consider the fiducial range ( 3.12 ) the 
numerical ef fect o f a(v) se ems to rem ain r ather small compared to the 'non ambiguous' contribution. Indeed, a{y) 
enters Eqs. ( 4.12 ), ( 4.14 ), ( [1.16] ), and ( 4.18 ) only through the differences i\(y) — cr(i>), ii(y) — <t(z/), and — a[y)\ 
i\{v) varies between 38.50 when v — 0, and 29.88 when v = 1/4, varies between 52.50 when v = 0, and 36.49 

when v = 1/4, while £5(1/) varies between 75.75 when v = 0, and 48.33 when v = 1/4. The a ambiguity therefore 
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has only a limited fractional effect on the coefficients it influences: < 10.0%, < 8.2%, and < 6.2% for the coefficients 
influenced through i\{v) — &(l>), %2{v) — <r(v), and i§{v) — o(y), respectively. 

By contrast we wish to emphasize that the use of the approximate ansatz (a la Wilson-Mathews |l(|) reducing 
the spatial metric to being conformally flat affects in a numerically much more drastic way the dynamical invariants. 
Indeed, the difference starts at the observationally much more significant 2PN order. The Wilson-Mathews truncation 
of the 2PN Hamiltonian is simply obtained by setting hJ T = in the formulas of JhJ . It reads (in ADM coordinates) 



ffw M (r,p) = E- - i + 1 {±( 3 „ - 1)(P 2 ) 2 - \ [(3 + ^)P 2 + Kn • P) 2 ] \ 



1 

2^2 



l(l-5, + 5, 2 ) (p»)S + jj(i-4*)(p») a ± 



4[(10+19,)p 2 -3Kn.p) 2 ]l-i(l + ,)l 



(4.22) 



Rieth [[L7| has computed (following, as above, the route of Ref. |14]]) the invariant functions of i?wM, Eq. ( 4.22 ). They 
read (see Eqs. (5) and (8) in 0) 



P 



1-kGM (-2E) 3 / 2 



11/ s 1 



13/ 1 

cr 2 I c 



~(60 - 181/ - 41i/ 2 ) ( 2E ) m - H( 7 - - - 2-V J ) E 1 



1 



1 



.16 j 
The corresponding Delaunay Hamiltonian reads 



— (140 - 54z/ - 31^W + — (60 - 18z/ - 41i/ 2 ) E 



24 



(4.23) 
(4.24) 



1 



jn 4 n z 



+ - 



; (140- 54i/— 31i/ 2 



j 3 n 



27 

j n 



- i(420- 42 l /-41i/ 2 )^ T 
8 V ; jn 3 



(29-3^-3^)-j 



(4.25) 



When comparing Eqs. flj.23| )-( [4.25|) to the 2PN-accurate versions of Eqs. ( |4.14|) , (|4.16|) , and (|4.18|) we see that 
the maximum fractional difference in the 2PN-level coefficients (between the exact -H2PN and the Wilson-Mathews 
truncation H WM ) is about 23% (for v = 1/4, and the coefficient of S 2 on the R.H.S. of Eq. ( 4.23| )). 2PN-level effects 
are crucial in determining important observables, such as the location of the Last Stable Orbit (as discussed, e.g., in 
Refs. p|,|T8|], and in our follow up work j^]), and are numerically more important than 3PN-level effects (even near 
the Last Stable Orbit the PN expansion parameter v 2 /c 2 ~ 1/6 is still smallish). Therefore we consider that the 
Wilson-Mathews truncation is, a priori, unacceptably inaccurate. W e sha ll come back below, and in Ref. to the 
more delicate problem of judging whether the 3PN-level uncertainty (3.12) on a is practically acceptable or not. 



V. DYNAMICAL INVARIANTS FOR THE 3PN CIRCULAR MOTION 



In the previous section we considered the invariant functions associated to generic (bound) 3PN orbits (E < 0). 
These orbits correspond to a doubly periodic rosette motion. The case of circular orbits, though quite special from a 
general point of view, plays a physically very important role and deserv es a dedicated study. Circular orbits represent 
the degenerate limit where the integration range [r m i n ,r max ] in Eq. (4.5) shrinks to a point. [In invariant terms, 
circular orbits are characterized by the fact that the motion becomes simply pe riod ic, instead of doubly periodic] 
Technically they are therfore characterized by setting the radial action i r , Eq. (4.5), to zero. In other words, the 
reduced 'principal quantum number' n = i r + j should be equalled to j. Using Eq. (4.18) this gives the following link, 
between energy and angular momentum, characterizing circular orbits: 



TP. . — 



1 

2f 



1 1 



1 



1 1 



c 2 4 



c 4 8~ 



v 



2 )-T + ^%M-^))^ 



(5.1) 
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where 



01 (v) 



3861 
128 



64 1 4l7r 



8833 \ 

— r 



313 
192 1 



128 



(5.2) 



Note that, evidently, it was not needed to derive the full Delaunay Hamiltonian ( 4.18]) to get the result (5A). It would 
have been sufficient to impose that the effective radial potential R(r,E,j), Eq. (4.4), have a double zero, i.e. that 
dR/dr = = d 2 R/dr 2 , or equivalently that the Hamiltonian H(r, p) (reduced to circular orbits by setting n • p = 0, 
p 2 = j 2 /r 2 ) has a minimum when considered as a function of r for fixed j. Let us also give ( though it is not an 
invariant function) the expression of the coordinate radius r (which is the shifted r' of Eq. (2.23)) in function of j 



1 4 



1 1 



c 2 j 2 c 4 8 y 'f c 6 



55 



163 2 1655 „ \ /61 „ 

-^-4 Wstatic j,+ (--5. k 



tic I 



An important observational quantity is the angular frequency of circular orbits, w c i rc . Let us fi rst recall that cj C i rc 
differs from the above considered periastron-to-periastron or radial frequency w ra dial> Eq. ( 4.20 ). In fact, because 
the frequency w ra diai mea sures the frequency in the 'rotating frame' of an orbit which precesses in space with the 
frequency Wpcriastron ( 4.21 ), it is easy to see that 



^radial ~t~ Wperiastron — 27T 



1 + k 
P ' 



Remembering that n = j along circular orbits we see from Eq. J4.2CD that 

1 dEr-.wc. 



GM dj 



(5.4) 



(5.5) 



which corresponds to the usual link d£ = oj c i TC dJ between the (unreduced) total energy £ — Mc 2 + £ NR and the 
total angular momentum J = [iGMj. As usual, it is convenient to introduce the (coordinate- invariant) dimcnsionlcss 
variable 



\ c 3 



2/3 



In the following, we shall set c = 1 to simplify the formulas. Inserting Eq. (5.1) into Eq. (|5.5|) gives the link 



GMlj c 



r 3/2 = ± 



1 + ^(9 + ^1 + 1(81-7^+ v 2 )^ + 8( 01 H - a{v))^ 
z J J J 



(5.6) 



(5.7) 



By inverting this relation one gets j as a function of x: 



7 ■ - x- 1/2 

Jcirc — ^ 



1 



1 



1 + -(9 + v)x + —(81 - hlv + v 2 )x 2 + -{o 2 {v) - <j{y))i 



where 



405 
128 



1 

64 



2 6889 \ 

4171^ - 1 



() 



r 



421 
192' 



3456 



(5. 



(5.9) 



At this point it is useful to consider the test-mass limit [y — > 0). Indeed, in this limit, i.e. in the case of a test 
particle m.2 moving, on a circular orbit, in a Schwarzschild background of mass mi (with mi <C mi), we have the 
relations (see, e.g., Eqs. (5.19)-(5.22) in pf) 



1 



J 



yj x(l — 2>X) 

x(l — 3x). 



(5.10a) 



(5.10b) 



Here the variable x is defined as in Eq. ( J5 . 6| ) , with M — mi + mi, and j = J/(fiGM) = J/(Gm\mi), where J is 
the total angular momentum of the system (which is entirely carried by the test particle). [At this stage we could 
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equivalently consid er that M is the mass of the black hole and [i = mim2/M t he m ass of the test particle.] The test 
mass result (|5.10h| ) suggests to focus on the function j 2 (x), which, from Eq. ( p.8|) , is given, when v ^ 0, by 



1, \ 25 2 16 , , N , , 
1 - -(9 + v)x + —vx 1 - —\pziy) - 



where 



03(f) 



1 

64 



41tH 



5269\ 



6 



- r 



511 
192 1 



1 

432* 



Let us also consider the energy for circular orbits. By inserting Eq. (5.11) into Eq. ( |5.1[ ) we get 

£total _ M 



E c 



1 

= X 

2 



1 - "^(9 + v)x - ^(81 - 57* + v 2 )x 2 - H ( 02 (u) - a(v))x 3 



(5.11) 



(5.12) 



(5.13) 



where 02(f) is given in Eq. ( |5.9| ). It is easily verified that Eqs. (5.8) and ( 5.13] ) satisfy (at the 3PN accuracy) the 
exact identity following from Eq. (|5.5|) 



dE, 



dx 



^3/2 djciic 

dx 



(5.14) 



When comparing Eq. ( 5.1 3| ) with the test-mass limit, there is a problem of definition of the best energy variable to 
use in a binary system. Indeed, if we consider a test mass of mass m 2 moving around a Schwarzschild black hole of 
mass mi , the total energy of the system reads 



£ 



total 



mi + £ 2 — Tax — k^p 2 — mi 



(5.15) 



where p 2 is the 4-momentum of the test mass m 2 , and £ 2 = —k^p 2 1S ^ ne conserved relativistic energy. Here, fc M is 
the time-tr ansla tion Killing vector which is (approximately) fc M = pi/mi, where p± is the 4-momentum of the large 
mass. Eq. (5.15), and the known results for circular motion in Schwarzschild (still with /1 = m\m 2 / {m\ + m 2 ) ~ m 2 
modulo 0{v)), yield 



E 



£ 



total 



„, „, cTM 

m i ~ m 2 ££ 



m 2 



rn 2 



1 - 2x 



(5.16) 



As emphasized in Ref. j| the expression ( ^5.15 ) is very asymmetric with respect to the labels 1 and 2. It seems much 
better to consider the Mandelstam invariant s = (£ tota1 ) 2 = —{p x + p 2 ) 2 = m\ + m 2 — 2p\ ■ p 2 and therefore the 
combination 



^£total^2 



pi ■ Pi £ 



TM 
2 



2m\m 2 



m\m 2 



m 2 



1 - 2x 
VI -3x' 



(5.17) 



This consideration, and the fact that the energy combination defined by the left-hand side of (5.17) still exhibits a 
branch cut singularity in the complex x-plane, motivated Ref. (J] to introduce a new invariant energy function e, 
defined by 



(gtotal)2 _ m 2 _ m 2 

2mim2 

In terms of E = £ NR = £ NR //i = (£ total — M)//x, the new energy function e reads 



(1 + E + -E' A y - 1 = 2E+ (1 + f)£ 2 + f£ 3 



£ 4 . 



Inserting Eq. ( pTl3| ) into Eq. ( |5T9| ) yields (at 3PN) 

e(x, v) = - 



1 - 3(3 + v)x - -^(36 - 35f)x 2 - y (04(f) - cr(f))x 3 



(5.18) 



(5.19) 



(5.20) 



13 



where 



, , 27 1 /41 , 4309 



77 

To 1 



270 



From Eq. (5.17) the test-mass limit of this function takes the simple (meromorphic) form 

1 - Ax 

e(ar,0) = 



1 - 3x 



(5.21) 



(5.22) 



Fo r co mpleteness, let us note that the test-mass limit of the invariant link between energy and angular momentum, 
Eq. (5.1), reads 



\/2 2 + JA(j) A , , 12 
l + E = ^ V W ; A(j) = l--. 



(5.23) 



Actually, the result ( 5.23| ) holds only along the sequence of stable circular orbits (of Schwarzschild radius R > 6GM, 
i.e. for a frequency parameter x = GM/R < 1/6). When x > 1/6 (but x < 1/3), i.e. for unstable circular orbits vA 
should be replaced by — VA in ( 5.23 ). This corresponds to the fact that the curve E — E(j) has a cusp at the Last 
Stable Orbit. 

To complete our list of exact results in the test-mass limit, let us also note that the periastron parameter 1 + k = 
<&/ (27r) is given, in the limit of circular orbits, by 0,^9) 

-1/4 

1 + k = I 1 - ^ ) . (5.24) 



12 

3 



Inserting Eq. ( |5.l|) into Eq. ( 4.24 ) we get, for circular orbits, 

3 



/■•-..< = -i + ^(45 - 12^1 + 6(05^) - 4, 
J * 3 3 



where 



, . 135 /ll 

OsM = h 7T^ !/ H 

5W 4 V 64 3 / 24 



101\ 53 2 
24 



From Eq. (|5.25|) follows that 



where 



(1 + fceire) -4 = 1 - j + 24^1 - 24(o 6 (^) - 



.'41 2 56\ 53 2 
bV y \64 3 / 24 



(5.25) 



(5.26) 



(5.27) 



(5.28) 



We will use in a further work S the various invariant functions computed above, at the 3PN level, to discuss the 
observable quantities associated to the Last Stable (circular) Orbit of a binary system. For the time being, we shall 
only remark that, contrary to what happened above in the dynamical invariants for generic orbits, one anticipates that 
the incomplete knowledge of the ambiguity parameter cr(y) might be muc h more serious in the dynamics of circular 
orbits. Indeed, if we consider for instance the function j~ 2 (x), Eq. J5.ll ), we see that u{v) modifies the coefficient 
of x 3 through the difference 03(f) — o"(f); 03(f) vanishes when v — and decre ases monotonically with v , taking the 
value —1.683 when v = 1/4. Therefore the addition of a(v), within the range (|3.12|), can modify a lot the coefficient 



of x 3 in j 2 (x). Similarly, in the case of the original energy function E{x), Eq. ( p. 13 ), <j(f) is subtracted from the 



coefficient 02(f) which decreases monotonically with v and varies between 02(0) = 3.164 and 02(1/4) = 0.397. This 
issue will be discussed in detail in 
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APPENDIX A: AMBIGUITIES IN THE 3PN ADM TWO-BODY POINT-MASS HAMILTONIAN 



In this appendix we recalculate the 3PN ADM two-b ody point-mass Hamiltonian given in Ref. 0] and we discuss 
the plausible range of the ambiguity parameter a, Eq. ( |3.8| ). [We use here the notation of 0), notably 16irG = 1.] 
The motivation for recalculating the Hamiltonian is three-fold: (i) We want to avoid using explicit solutions of some 
Poisson equations with complicated (and singular) source terms (it was shown in Ref. [|l3| that it may not be possible 
to obtain unique solutions to these equations), (ii) Before applying any regularization procedure we want to represent 
the Hamiltonian as a volume integral with an integrand which diverges locally as weakly as possible (thus absorbing 
stronger divergent terms into surface integrals), (iii) We wish to eliminate all 'contact' terms (Dirac delta functions) 
in the Hamiltonian, so that only one regularization procedure is needed (the one used to regularize the volume integral 
giving the Hamiltonian). 

Our starting point is the 3PN Hamiltonian H 3 p^ as given by Eqs. (54)-(60) of Ref. In the first step we eliminate 
from H 31 (given by Eq. (55) of j|) the function (8)1 by integrating by partsQ (8)1 ^ Q m <5 Q = -0(8)1 ( A 0(2)) = 
— (A0(8)i) 0(2)+ exact divergence. To do this we use the equation (see Eqs. (All) and (A12) in 

fl 1 1 f 5 [j; 9 \ p 2 

A0( 8 )i = ^ { FTo^(2) " T7770(2)0(4) + J0(6) + 7770(4) ~ 7777 



t \512^ 2 ) 32^,^4, Tgvw-r ^^(4, 12g n 2 

jLj, (Pa) 2 1 (Pa) 3 1 PaiPaj h TT 

64 9{2) mi 16 m« 2 ml {i)rj 



| m a 6 a + 7^0(4)^^ + \ A {( h Wv) 2 } ■ ( A1 ) 



When performing an integration by parts we always check that the exact-divergence terms fall off at infinity faster 
than 1/r 4 , so that (formally) they do not contribute to the Hamiltonian. In the next step we split the function 0( 6 ) 
into three parts (see Eq. (A4) in ||]) 

0(6) = 0(6)1 + 0(6)2 + 0(6)3, (A2) 

and we eliminate the functions 0(6)2 and 0(6)3- To perform this we use again integration by parts and equations (see 
Eqs. (A6) and (A7) in §) 

A0 ( 6)2 = -(^ 3 )) 2 > (A3) 



1 



A 0(6)3 = 770(2),y/l(4)y' (A4) 

The elimination of the functions 0(g)i, 0(6)2, an d 0(6)3 means that the goal indicated in the item (i) at the beginning 
of this appendix is achieved; the Hamiltonian depends only on the three lowest-order solutions of the Hamiltonian 
constraint: 0(2), 0(4), and 0(6)i- These three functions are not influenced by the possible ambiguities in solving the 
Hamiltonian constraint discussed in Ref. Jl3| . To achieve the goals of items (ii) and (iii) one must perform many 
integration by parts and employ several times the equations fulfilled by the functions 0(2), 0(4), and 0(6) i, which read 
(see Eqs. (15), (16), and (A5) in §) 

A0( 2 ) = - ^m a S a , (A5) 

a 

%)=Efo(2)-^)«A. (A6) 

a \ a / 

A 0(6)i = E (-h?<a + + 7^( 2 )S + 5%r) maS - ( A7 ) 
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a ~ ~ . 

Finally, we are able to put the 3PN two-body point-mass ADM Hamiltonian in the following form (dropping many 
surface terms, after checking that they fall off fast enough at infinity not to contribute to the Hamiltonian) 



7 In this appendix we abbreviate 5(x — x a ) by Sa- 
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H^PN — — - 
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/ d 3 x(h 1 + h 2 + h 3 ), 

n a •> 



(A8) 



where 



kl ~ 8 {~k^ {2) ~ T^ (2)<?!>(4) + <?!>(6)1 ) S[A) > 1 + 16 (~8~^ 2) ~ 3< ^ (4) ) 5(4) ' 4 + f^ 2 )'^ 8 )'* ~~ i^Wy.fe^ 6 )^'^ 
+ G 5(4) ~ 4^ (4) ) fe) + (iV (4)j ~ I^W-j) ^( 2 ).* ft Wti ~ ((^(2)^(3)) 
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In Eqs. ( A9a )-( |A9c|) the following notation was introduced 



5, 
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(A9c) 

(AlOa) 
(AlOb) 
(AlOc) 
(AlOd) 
(AlOe) 



We have not found any ambiguity connected with the regularization of the divergent integrals contained in the first 
part hi of the Hamiltonian (see below for details). Regularization of the hi part contributes only to the static 
ambiguity, and regularization of the fa part contributes only to the kinetic ambiguity. The static part of the 3PN 
Hamiltonian (discussed in detail in Ref. |D|; see footnote |^ in the present paper) is equal to j d 3 x h 2 - 

We regularize all the divergent integrals involved in the Hamiltonian given by Eqs. (|A^)~( |A9[) by means of the 
Riesz-formula-based regularization technique described in Appendix B 2 of j| and rediscussed in Sec. Ill of the 
present paper. The final result, after elim inat ing the time deri vativ es of the particle positions and momenta by using 
the Newtonian equations of motion, Eq. ( |3.l| ), is given in Eq. (3.6). 

To investigate the lack of u niqueness of the regularization method we have used the following procedure. To many 
individual terms in Eqs. (A9) one can associate exact divergences which must be adde d to these terms to replace 
them by some of their integration-by-parts equivalents. E.g., with the first term in Eq. (A9a) one can associate two 
exact divergences, — y|g [(0( 2 )).i<5(4)].i an d — i§g [0(2)^(4), The result of the regularization method is stable against 
integration by parts only if the regularized value of the integral of any of these exact divergences is zero. 
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With hi, Eq. (A9a), one can associate the following 18 exact divergences 
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1 



(6)l5( 4 ),i] i , 
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(All) 



The question now is whether, when explicitly computed (and regularized) by the method we are going to recall, the 
integrals of the exact divergences in the list ( All| ) vanish or not. The explicit method of regularized- integration we 
use is the following: First, we explicitly perform the differentiations present in the divergences using, if necessary, the 
distributional rule of differentiation of homogeneous functions described in Appendix B 4 of 0] (cf. the example given 
there). Note that we always use Leibniz's rule when differentiating products of functions, the distr ibut ional rule is 
used only for individual functions (by which we understand the functions explicitly enter ing Eqs. (A9): 0(2), 0(4), 



(6)1, 



?(4), . . . ; it implies that in computing the new reference Hamiltonian given by Eqs. ( |A9| ) we have not used the 
distributional rule). E.g., we compute (0(2)0(2),j) m two steps: (i) using Leibniz's rule one gets O?2)0( 2 ),j 



20(2) 0(2) ,i 0(2) ,j + 0(2)0(2), y i (ii) the distributional rule is applied only to (/){2),ij- After this, a typical integral consists 
of two parts: without and with Dirac deltas. The first part is computed by means of the Riesz-formula-based 
regularization technique, while, for the second part, we have used the Hadamard-partie-finie regularization described 
in Appendix B 1 of g|. In the second case we have also investigated the stability of the result against 'threading' the 
partie finie (Pf) over a product of functions, i.e., we have checked whether 

Pf(/ 5 ) = Pf(/)Pf( 5 ). (A12) 

For all divergences in (All) the final result of these explicit computations is zero (regardless which side of Eq. ( A12| ) 
is employed in the contact terms). 

With h,2, Eq. (A9b), one can similarly associate the following six integrals of exact divergences 
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(0(2)) .0(6)1 



^(2)0(4)), 0(4) 
d 3 X [0(2) 0(4) 0(4), i] ; 



(f.l 



(A13b) 
(A13c) 
(A13d) 
(A13e) 



A 26 = - 



(fx 



(Al3f) 



Below we denote by primes the result of regularizing with the use of the left-hand side of Eq. ( A12 ) (for these parts 
of the integrands whic h are proportional to Dirac deltas) , and by double primes the result obta ined when using the 
right-hand side of Eq. (A12). With this notation, the regularized values of the integrals in (A.13) read 
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With Eq. ( A9c ), we associate eight integrals of exact divergences 
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(A14a) 

(Al4b) 
(Al4c) 

(Al5a) 
(A15b) 
(Al5c) 
(A15d) 
(Al5e) 
(Al5f) 
(Al5g) 
(A15h) 



The results of the regularized-integration of ( A15| ) read (here k = v 2 (p 2 — 3(n • p) 2 ) /r 3 ; note that the exact diver- 
gences A32 and A34 do not contain contact terms) 
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(Al6a) 
(A16b) 
(A16c) 
(A16d) 
(Al6e) 

(Al6f) 

(Al6g) 
(A16h) 



From Eqs. (A14) and (Al£) follows that only one out of the 32 exact divergences considered above, A21 = vj(Ar ), 
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is connected to the static ambiguity parameter ^static, whereas eight exact divergences A31-A38 'contribute' to the 
kinetic ambiguity parameter akinetic; their regularized values are between —10.9k and +8k. 

We are not confident that the so- found relatively small contribution (±1/4) to w s tatic is indicative of a smaller 
ambiguity in cj sta tic than in Wkinotic- Therefore, we used also another (more primitive, but hopefully more robust) way 
of estimating the plausible values of the ambiguity parameters e cstat ic and ^kinetic- We simply looked at the values of 
all the numerical coefficients in the final 3PN Hamiltonian, Eq. ( |3.6|) . The Hamiltonian has 27 numerical coefficients, 
entering as factors of the different monomials made of powers of p 2 , (n • p) 2 , 1/r, and v. Out of them only three 
are influenced by ambiguities (the coefficients of v 2 p 2 /r 3 , v 2 (x\ ■ p) 2 /r 3 , and v/r A ). The numerical values of the 
unambiguous coefficients range between 7r 2 /64 — 335/48 w —6.825 (i/p 2 /r 3 ) and 17/2 = 8.5 (y (p 2 ) 2 /?" 2 ), whereas 
the ambiguous coefficients in the reference Hamiltonian (fT(j) are equal to —55/12 w —4.583 (i^ 2 p 2 /r 3 ), 27/8 = 3.375 
(i/ 2 (n • p) 2 /r 3 ), and 109/12 - 21tt 2 /32 w 2.606 {v/r 4 ). 

As a result of the above discussion we take, as plausible ranges for the ambiguity parameters astatic and akinetic, 
the ranges 

- 10 < Static < 10, -10 < kinetic < 10. (A17) 



APPENDIX B: COMPUTATION OF THE RADIAL ACTION INTEGRAL 



The explicit form of the coefficients A, B, C, Di entering the effective radial potential R(r,E,j), Eq. (4.4), reads 
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(Ble) 
(Blf) 
(Big) 
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(Blh) 



The PN expansion of the effective radial action integral i r , Eq. (4.5), can be writt en (as explained in Sec. IV of 
the present paper) as a linear combination of the integrals J q , m defined in Eq. (4.1C). The only integrals that one 
needs to compute at 3PN order are J q . m with < q < 3: when q — 0, to — 0; when q = 1, 3 < m < 7; when q — 2, 
6 < to < 8; and when q = 3, to = 9. The explicit results of the integration (performed by the two methods indicated 
in the text) read 
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^1,3 = — ^TT, (B2b) 



(-C)i 
IB 2 - 

5B 2 - 

7(-c)i 

70AB 

2(-C)\ 

-7B 2 - 
2(-C)i 

+ 

"8(-rr 

2(-Cf 



35 2 - AC 

Ji 4 = =— , (B2c) 

2(-C)i 

Ji,5 - 2( _ C)J , (B2d) 
35B 4 - 30AB 2 C + 3A 2 C 2 

J " " iP^ • (B2c) 

S(63B 4 - 70AB 2 C + \hA 2 C 2 ) 

■A,7 = — — n 1 (.B21J 



_ 3(-5i3 2 + AC) 

J2 ' 6 - 2(-C)i ' ( g) 

_ 5£(-7g 2 + 3.4C) 

J2 ^ 2^)1 ' (B2h) 

15(-21B 4 + 14AB 2 C* - A 2 C 2 ) . , 

^,8 = > 

^= 35 ^ C) (B2j) 



[1] T. Damour and N. Deruelle, Phys. Lett. 87A, 81 (1981); T. Damour, C. R. Seances Acad. Sci., Ser. 2 294, 1355 (1982). 
[2] P. Jaranowski and G. Schafer, Phys. Rev. D 57, 7274 (1998). 

[3] L. Blanchet, T. Damour, B. R. Iyer, C. M. Will, and A. G. Wiseman, Phys. Rev. Lett. 74, 3515 (1995); L. Blanchet, T. 
Damour, and B. R. Iyer, Phys. Rev. D 51, 5360 (1995); C. M. Will and A. G. Wiseman, Phys. Rev. D 54, 4813 (1996); L. 
Blanchet, B. R. Iyer, C. M. Will, and A. G. Wiseman, Class. Quantum Grav. 13, 575 (1996); L. Blanchet, Phys. Rev. D 
54, 1417 (1996). 

[4] T. Damour, B. R. Iyer, and B. S. Sathyaprakash, Phys. Rev. D 57, 885 (1998); (work in progress). 
[5] T. Damour, P. Jaranowski, and G. Schafer (in preparation). 

[6] L. Blanchet and G. Fayc (work in progress); M. E. Pati and CM. Will (work in progress). 
[7] T. Damour and G. Schafer, Gen. Rel. Grav. 17, 879 (1985). 
[8] G. Schafer, Phys. Lett. 100A, 128 (1984). 

[9] T. Ohta, H. Okamura, T. Kimura, and K. Hiida, Prog. Theor. Phys. 51, 1220 (1974). 
[10] G. Schafer, Ann. of Phys. (N.Y.) 161, 81 (1985). 
[11] B. M. Barker, and R. F. O'Connell, Can. J. Phys. 58, 1659 (1980). 

[12] T. Damour, in Gravitational Radiation, edited by N. Deruelle and T. Piran (North-Holland, Amsterdam, 1983), pp. 59-144. 
[13] P. Jaranowski and G. Schafer, Phys. Rev. D 60, 124003 (1999). 
[14] T. Damour and G. Schafer, Nuovo Cimento B 101, 127 (1988). 

[15] A. Sommerfeld, Atombau und Spektrallinien, 8th ed. (Vieweg, Braunschweig, 1969), Vol. 1. 

[16] J. R. Wilson and G. J. Mathews, Phys. Rev. Lett. 75, 4161 (1995); J. R. Wilson, G. J. Mathews, and P. Marronetti, Phys. 
Rev. D 54, 1317 (1996). 



20 



[17] R. Rieth, in Mathematics of Gravitation. Part II. Gravitational Wave Detection, edited by A. Krolak, Banach Center 

Publications (Institute of Mathematics, Polish Academy of Sciences, Warszawa, 1997), Vol. 41, Part II, pp. 71-74. 
[18] A. Buonanno and T. Damour, Phys. Rev. D 59, 084006 (1999). 
[19] W. Junker and G. Schafer, Mon. Not. R. Astron. Soc. 254, 146 (1992). 



21 



